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The objective of this research is to study small transition metal clusters 

at a high level of approximation i.e. including all the valence electrons in 

the calculation and also including extensive electron correlation. A major 

goal of this study is to understand the electronic structure of these small 

metal clusters and by comparison of dimers, trimers and possibly higher 

clusters to extrapolate the information obtained to provide insights into 

the electronic structure of bulk transition metals. Small metal clusters 

are currently of considerable experimental interest and some information is 

currently becomming available both from matrix ESR studies and from gas 

phase spectroscopy. Collaboration between theorists and experimentalists 

Is thus expected to be especially profitable at this time since there is some 

> 

experimental information which can serve to guide the. theoretical work. 
It should be pointed out here that transition metal ifiolecules have been 
difficult for both theory and experiment and a close collaboration between 
theory and experiment is likely to lead to a deeper understanding than 
would be obtained by either approach alone. 

These studies are expected to be applied iu several areas. One area 
of current interest is the reactivity of small metal clusters. As an example 
Smalley et. al. have studied the reactivity of H2 with Nb clusters and 
find that Nb atom and Nb 2 are nnreactive to H 2 while Nb 3 is reactive 
to H2. Similarly, several groups have observed that transition metals form 
compounds of specific stochiometry with O or H atoms. These observations 
are relevant to understanding the catalytic properties of small clusters but 
the thoretical basis for understanding these observations does not exist 
because of the very limited understanding of the electronic structure of 
these metal clusters. In this area theory may well be able to provide results 


more rapidly than experiment. For example for the Nb 2 molecule theory 
is able to study the potential curves as easily as for the V 2 molecule but 
the spectroscopic analysis for Nb 2 is much more difficult than for V 2 and 
has not been possible to date. Similarly it should be possible to study Nb 3 
computationally and a comparison of the electronic structure of Nb 2 and 
Nb 3 should provide some insight into why one is reactive to H 2 while the 
other is not. 

A second area of application of these studies is in understanding the 

interaction of bulk metals with hydrogen i.e. hydrogen embrittlement. 

Here studies had been carried out earlier of the potential function for a 

single hydrogen atom interior to BCC iron using a one-electron ECP which 
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replaces the 3d 7 Fe atom configuration as well as th&.Ar core Le. only 
a single 4s electron is left on each iron atom (Which* Work carried out 
under NCC2-148). Here the plan is to place a small cluster of iron atoms 
with the 3d electrons included inside a larger cluster with the remaining 
atoms treated at the one-electron ECP level. This approach would permit 
studying features which would not be treated reliably by the one-electron 
ECP approach such as the formation of an FeH 2 phase where metal 3d 
electrons may be involved in the bonding. 

in the first six months of this grant calculations have been carried out 
for the the Y 2 , Nt> 2 , and Mo 2 molecules thus extending the earlier studies 
of the first row transition metal dimers to elements of the second transition 
row. These studies show that the second transition row elements form 
stronger nd bonds due to i) a larger 4d to 5s orbital stye ratio as compared 
to the 3d to 4s orbital size ratio for the first transition row and ii) smaller 
nd exchange interactions for the second transition row as compared to the 


first transition row. Calculations were also carried out for the TiV molecule 
which was of experimental interest since it has been observed via ESR in a 
matrix. Here the calculations confirmed the 4 £"" ground state symmetry 
which had been deduced from the ESR experiments. The details of this 
work along with the previous work on first row transition metal dimers have 
been presented in a paper which is included as an appendix to this report 
and will not be discussed further here. 

In addition to the studies of diatomic molecules presented above studies 
are currently in progress for selected transition metal trimers. Currently 
the Ca 3 , Sc 3 , Y 3 , and Cu 3 molecules are being studied. - • 

For the Cu 3 molecule all-electron studies are beingvcarried out for the 

2 E state arising from three 4s 1 3d 10 Cu atoms. These studies are single 

% 

configuration SCF plus configuration interaction (Cl) calculations with an 
[8s6p3d] basis set. From studies on CU 2 one expects these calculations to 
lead to somewhat too long Cu-Cu bond lengths due to a combination of 
relativistic and size consistency effects. However, these calculations should 
lead to accurate relative energetics for the various geometric isomers of Cu 3 . 
Two magor structural possibilities exist here, linear and near equilateral 
triangle. Because the equilateral triangle form of Cu 3 is a 2 E state one 
expects a Jahn-Teller distortion leading to 2 Ai and 2 B 2 states in C 2v 
symmetry. For the near equilateral triangle geometries the Cu-Cu distance 
was first optimized for equilateral triangle geometries leading to an optimal 
Cu-Cu distance of 4.60 &o. The 2 Ai and 2 B 2 states were then studied as a 
function of Cu-Cu-Cu angle with two sides fixed at 4.60 ao. The resulting 
curves are shown in Fig. 1 where one sees that the 2 E state is a maximum 
on the potential surface while the 2 Aj state has an aclite angle geometry 





and the 2 B 2 state has an obtuse angle geometry. The details of the potential 
surface interconnecting these isomers are under study. One possibility is 
that the 2 B 2 state is one of three equivalent minima while the 2 Ai state is 
one of three equivalent saddle points. The calculations do indicate that 
the linear isomer of C 113 is at least 0.2 eV higher than the equilateral 
triangle geometries. This result is consistent with the interpretation of 
the spectroscopic results of Smalley et. al. where it is concluded that the 
ground state of C 113 is a strongly Jahn- Teller distorted molecule. 

Calculations are also being carried out for SC3 and Y3. These systems 
were selected for study because of the ESR studies of Weltner et.al. which 
establish that the ground state of Sc 3 is a doublet state with equivalent 
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Sc atoms while Y 3 is a doublet state with inequival^nt Y atoms. The 
calculations are made tractable by the use of an ECP tcf "replace the Ne core 
for Sc and the Ar core for Y. This leaves the ns and np core levels in the 
valence space and is more reliable than had these levels been incorporated 
into the ECP. The ECPs were developed by Wadt and Hay and have been 
tested against all-electron calculations for Sc 2 , Y 2 , and Nb 2 . 

The calculations are CASSCF calculations with the (n-j-l)s and nd 
in the active space and use a [4s4p3dlf] valence basis set. With this size 
active space for Sc 3 there are a small enough number of electrons that 
we can allow the MCSCF procedure to scan all possible states. The first 
encouraging result is that the doublet states are lower than the high spin 
states and the wavefunctions show very definite indications of 3d bonding. 
Three doublet states have been studied. The dominant configurations for 
these states are: 


2 A' 4suaf 4 sue ' 4 3d Sd^a '/ 2 
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2 A" 4saa' 1 2 4 sae' 4 3da&f 

2 E 4saa , 1 2 4sae' 3 3dad!£ 3d7ra'/ 2 

Figure 2 shows the energy of the 2 A', the 2 A" and the 2 E states as a 
function of the Sc-Sc distance for equilateral triangle geometries, The 2 E 
state dissociates to two ground state atoms and one excited state atom while 
the other two states dissociate to three ground state atoms. The CASSCF 
curve for 2 E shows a well near R Sc-Sc of 5.5 au, while the 2 A" shows only 
repulsive behavior as far out as 5.75 au. Preliminary contracted Cl results 
indicate that the 2 A' and 2 A" states are bound by about 0.75 eV at R Sc-Sc 
of about 5.5 ao and also one expects the 2 E state to be bound at a geometry 
near the CASSCF geometry. However,the corresponding contracted Cl 
calculation for this state may be beyond present computational capabilities. 

K 

The state that seems most in line with the ESR results is the 2 A' or 2 A" 
state since these states would not Jahn-Teiler distort and have almost no 4s 
character in the singly occupied orbital. Also there is a very low-lying state 
observed by Moskovitts et.al. which is consistent with the small separation 
between the 2 A' and 2 A" states. 

If the true ground state is 2 A' or 2 A" that implies some real differences 
between Sc 2 and Sc 3 since Sc 2 has a ground state which arises by 4s 
to 3d promotion while the 2 A' state is bonding largely because of 4s to 
4p promotion. The populations here are 4s=4.50, 4p=1.21,3d=3.25,and 
4f=0.04 for the 2 A' state at R Sc-Sc of 5.5 ao. 

Calculations have also been carried out for Y3 at the CASSCF level. 
Here the lowest configuration at the CASSCF level for. equilateral triangle 
geometries is the 2 A' state.(See Figure 3.) This state would not be con- 
sistent with the nonequivalent Y atoms observed in the ESR spectrum. 
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However, for linear geometries there is a low-lying 2 £+ state which could 
be consistent with the non-equivalent Y atoms, The populations for Y 3 are 
5s^3.69,5p^2.25,4d^2.99, and 4f=0.06 for the 2 A' state at R Y-Y-6.3 
ao which is the minimum energy for equilateral triangle geometries. The 
much larger valence p population here is consistent with a lower s to p 
promotion energy for the Y atom as compared to the Sc atom and may be 
consistent with a significant change of geometry between SC 3 and Y 3 . This 
point requires further study. 
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Appendix A, The following manuscript describes the calculations for transi- 
tion metal dimers, 
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Abstract 

The results of CASSCF calculations are presented for the SC 2 , Ti 2 , 

r 

V 2 , Cra, Cua, TIY, Ya, Nt> 2 , and M 02 molecules. CASSCF/CJ calculations 
were also carried out for Sc 2 , Ti 2 > Cu 2) and Y 2 . The CASSCF procedure 
is found to generally provide reliable R, and <*/« values, flowever, D # values 
are systematically underestimated at the CASSCF level;the worst case is 
Cr 2 where the CASSCF curve is not bound although a shoulder is observed- 
in the region near the experimental R e . The CASSCF procedure is shown 

to provide a consistent set of calculations for these molecules from which 

» 

trends and a simple qualitative picture of the electronic structure may be 
derived. These calculators confirm the 5 £~ ground state of Sc 2 and the 
4 E“ ground state of TiV and lead to predictions for other molecules in this 
aeries. So far only the 8 E“ ground state symmetry and the bond length of 
V 2 have been confirmed by experiment. 


fMailing address: NASA Ames Research Center, Moffett Field, CA 94035 


I, Introduction 


Transition metals and transition metaJ(TM) compounds are currently 

of considerable interest because of their relevance to catalysis and to materials 

* 

* 

science problems such as hydrogen embrittlement and crack propagation in 
metals. In spite of such interest, progress in understanding the chemistry of 
TM molecules has been slow due to both experimental and theoretical 
difficulties. For a review of the current state of both theory and experiment 
for TM molecules the reader is referred to the review article by Weltner and 
Van Zee[l]. Recent progress in the experimental characterization of TM 
molecules has come from both matrix isolation and gas phase spectroscopic 

approaches, while progress in the theoretical treatment of TM molecules has 
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occurred as a result of improved methods of treating electron correlation. 

These improvements in part involve more efficient computational methods 

* 

m 

such as improved MCSCF[2], CASSCF[3], and direct CI[4] techniques but 

* 

also qualitative ideas such ms the GVB[5] method have provided .a con- 
ceptual framework for understanding electron correlation effects in these 
molecules. 

A dramatic example of the importance of electron correlation for transi- 
tion metal compounds is given by the NiCO X X E+ state. This state is un- 
bound by 2.6eV at the SCF level, while the triplet states are only slightly 
repulsive [6].. The Mulliken populations for the triplet states show a 3d 
population of 9, while the , £+ state has a 3d population near 10. If an 
MCSCF approach is used, the 3d population for the triplet states is only 
slightly changed, but the 3d population for the *£+ state is qualitatively 
different, 9.5. This arises because of the mixing of the Ni •s'Sd’' (*D) and 
3d 10 (*S) states, and the l £+ state is now the ground state of NiCO. 


It would take a very large Cl to over come the 3d 10 bias of the SCF 
reference. In this case, the CASSCF wavefunetion yields results in remark- 
able agreement with a multi-reference singles and doubles Cl calculation [6] 
and a D e of 1.3eV which is quite consistent with experiment (1.3±0.6 eV 
[ 7 ]). 


MCSCF treatments do not always yield results as good as those ob- 
tained for NiCO but they usually do represent a m^jor improvement over 
SCJ results. For an especially difficult case, the Cr 2 molecule, SCF leads 
to a potential curve which is unbound with respect to the SCF atoms by 10 
eV in the region near the experimental R f . An MCSCF treatment shows 
the outer well arising from the 4$-4s bonding but Is still unbound by 1.4 eV 

in the region near R, where one expects multiple 3d-3d'bonds [8,9J. While 
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this might appear to be a major failing of the MCSCF approach, if Cr 2 is 
compared to Moa [8] and V 2 (9], both of which are bound at the MCSCF* 
level, useful information about the nature of the bonding in Cri can be in- 
ferred. Thus, unlike the chemistry of the first two rows of the periodic table, 
where computational chemistry is achieving accuracy comparable with ex- 
periment, for some transition metal compounds only qualitative accuracy 
is pcssibleJIowever, so little is known about many of these systems that 
even qualitative information can contribute to a basic understanding of the 
chemistry. This is somewhat pessimistic since the MCSCF procedure does 
lead to good R e and values for many of the TM dimers even if the bind- 
ing energies are systematically underestimated and can treat states arising 
from the same occupation to very similar accuracy thus in these cases T c 
values should be good. When comparing very different states, sometimes a 
combination of theory and experiment can combine to* resolve a question. 
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l For example for Ti 2 theory can calculate the R e and w, of the low-lying 

and 7 £+ states but can not clearly establish which iff the ground state. 
However, combining the calculated results with the experimental u e leads 
to a clear assignment of a 'Lf ground state and provides new information 
about the chemistry. 

Many people have investigated individual transition metal dimers (see 
review article of Weltner and Van Zee [1]). However, as noted for Cr 2 , 
a consistent set of calculations can reveal trends which greatly contribute 
to the understanding. Goodgame and Goddard have considered Cr 2 and 
Moa |8], Shim and co-workers [10-11] have considered several transition 
metal dimers for those metal atoms with more than half filled nd shells, 
and these authors have considered many of the dimers for those metal 
atoms with less than half filled nd shells where strong nd-nd bonding 
is common. In this article we summarize our previous work, report on 
several unpublished systems and describe the trends for the bonding in the 
transition metal dimers on the left side of the first two transition metal rows. 
Here we concentrate on the CASSCF results since this level of calculation is 
possible for all the systems and provides a useful qualitative understanding. 
However, for some systems we also are able to carry out multi-reference 
singles and doubles CI(MRSDCI) calculations using the CASSCF orbitals 
and these calculations serve to calibrate the CASSCF approach. For other 
systems, such as Cra, the CASSCF -j-MRSDCI approach is not possible. 
This has led to several attempts to understand the bonding using even 
more qualitative approaches [12-14]. These include the extended MCSCF 
approaches [12-13] and approaches based upon physical arguments [14]. 

Section II describes some trends for the transition .metal atoms which 
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are important for understanding the bonding of the TM dimers and also 
discusses qualitative features of the bonding in the TM dimers. Section 
HI discusses the technical details of the calculations. Section IV discusses 
the potential curves and spectroscopic constants obtained for the first row 
TM dimers 'while Section V describes the same information for the second 
row TM dimers and discusses some of the differences between the first 
and second row TM dimers. Section VI discusses the electron correlation 
problem for O 2 and discusses some possible future approaches that might 
be used. These future approaches are discussed in relation to some of the 
current more approximate techniques. Finally Section VII presents the 
conclusions from this work. 

IL Qualitative Features of the Bonding in the TM Dimers 

% 

We first consider some features of the TM atoms, in particular the 
relative (n+l)s and nd orbital sizes, the nd-nd exchange terms, and the or- 
dering and separations of the low-lying atomic states[15j. These properties 
of the atoms are important because they control the atomic states available 
for bonding, the strength of the nd-nd bonds and the nature of the bond- 
ing i.e. whether the bonding orbitals involve predominately (n+l)s or nd 
character. Table I shows the following properties i)<r( n+ i) e >/<r n d>, 
obtained form numerical Hartree-Foc’ calculations for the (n-f-l)s 1 nd rn+1 
state ii) the (n+l)s 2 nd m ->(n-bl)s 1 nd m "*' 1 excitation energy and iii) the 
(n4-l)s 2 nd m --^nd r7, + 2 excitation energy. Looking first at the relative orbi- 
tal sizes we see tbat for the first transition row the ratio <r 4 ,>/<r 3< j> 
increases monotonically from 2.03 for Sc to 3.36 for Cu. The same trend 
is seen for the second transition row but here the 5s and 4d orbital sizes 
are more comparable. The increase in <r(*+i),>/<r* < f> as one moves 
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across a row is due to increased shielding of the nuclear charge by the nd 
electrons. 

Looking next at the excitation energies, the ground state for the first 
few elements of the first and second transition rows is the (n-f l)s 2 nd m state, 

The trends in the (n-f l)s 2 nd m -*(n-H)s 1 nd m + 1 excitation energies depend 
on two competing factors i) the stabilization of the nd level with respect 

* 

to the (n-f l)s level as the nuclear charge increases and ii) the number of 
nd-nd exchange terms. XJp to the point that the nd shell is half filled both 
of these factors favor the (n-H)s l nd m + 1 state for larger nuclear charge 
since movement of an s electron into the d shell increases the number of 

nd-nd exchange terms. This leads to a monotonic decrease in the excitation 

V- 

energies from Sc to Cr and from Y to Mo. Starting with Mn or.-Tc, the 
excited state has one more doubly occupied nd orbital in each case than the 
ground state. Thus, s to d promotion leads to loss of nd-nd exchange terms 
and a strong preference for the (n-f l)s 2 nd 5 state resulting in a discontinuity 
between Cr-Mn and Mo-Tc. However, from Mn to Cu and Tc to Ag, there 
is once again a monotonic decrease in the excitation energy due to the 
stabilization of nd with respect to (n-f l)s for increasing nuclear charge. 

From Table I one also sees that the nd m + 2 state is an excited state for the 
first transition row but becomes a low-lying state for the right half of the 
second transition row. In fact 4d 10 is the ground state configuration for 
Pd. 

For the molecules considered here the atomic states whfcN are sufficiently 
low-lying to be available for bonding are the (n-f l^nd"* and (n-f ljs 1 nd m+1 
states. These two states lead to three different atomic asymptotes: i) 
(n-fl^nd" 1 -f (n-fl)s 2 nd m , ii) (n-f l)s 2 nd m -f (n-fl)s ! nd m + 1 and iii) 
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* (n-f"l)s , nd m+1 -f (n-H)s 1 nd m +*. The accessibility of these asymptotes 
of course depends on the (n-H)s 2 nd m -*(n+l)s 1 nd m + 1 excitation energy. 
Thus for the first transition row we find that for Sc 2 the low-lying states 
arise from the first two asymptotes with 4s , 3d m + 1 -f 4s , 3d m + 1 too high 
In energy to lead to the ground state, while for Ti 2 , V 2 , Cr 2 , and Cu 2 
the ground states arise from the 4s , 3d m + 1 ■+• 4s , 3d m+l atomic asymptote. 

For the first transition row, because the 4s orbital is significantly 
larger than the 3d orbital the predominant interaction in the TM dimers 
at large internuclear separation (R) is between the 4s orbitals with very 
little 3d interaction. For states arising from the 4s 2 3d m -f4s 2 3d fn atomic 
asymptote this interaction is basically repulsive and only a shallow well at 
large R f (approximately. 8.0 ao) arising from the 4s-*>4p near degeneracy 
effect is observed. For states arising from the 4s 2 3d f ” +4s , 3d m + 1 atomic 
asymptote the 4s interaction is weakly bonding at intermediate R but is. 
repulsive at small R leading to intermediate R e v»lues(approximately 5.0 
ao). At these R e values the 3d-3d overlaps(S) are small which favors one- 
electron 3d bonds, (ie. one electron per bonding orbital as in H^~, a bond 
order of one half) whose bonding terms vary with distance like S, over two- 
electron 3d bonds (ie. two electrons per bonding orbital as in H 2 , a bond 
order of one), whose bonding terms vary with distance like S 2 [17]. Finally 
for states arising from the 4s 1 3d m + 1 -f* 4s 1 3d m+1 atomic asymptote the 
4s~4s portion of the interaction is attractive and also appears to not become 
repulsive until well inside the optimal 4s-4s bonding radius. Thus, states 
arising from this atomic asymptote favor short R c regions where the 3d-3d 
overlaps are large enough to permit two-electron 3d bonding. 

The strength of multiple one-electron nd bonds is 'determined mainly 
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* by overlap considerations, since the resulting singly occupied orbitals are 
orthogonal and high spin coupled, leading to no loss of atomic exchange 
interactions. The relative overlaps of nd orbitals are nda p^ndOnd#. 
Thus, the nda and ndir orbitals are filled first. An example of such bonding 
is the state of SC 2 , which has three one-electron 3d bonds of a and n 
symmetry. 

The strength of multiple two-electron nd bonds is determined both by 
the overlap of the orbitals and by the size of the atomic exchange terms 
which are lost due to bond formation. Overlap considerations primarily 
determine the occupation of the bonding orbitals. Since two electron bonds 

vary with R as S 2 , and one-electron bonds as S, the nd# bonding is even less 

* 

a 

favorable relative to ndo and ndrr than for one-electron, .bonds. Thus, for 
Tfe the state which arises from 4s ! 3d a + ds^ 8 ha? three two-electron 
3d bonds with bonding orbitals of a and t symmetry. It is interesting that 
in the region where the ndo and nd* orbitals have large overlaps leading 
to two-electron bonds, the overlaps of the nd# orbitals are so small that 
one-electron bonds are favored. For example in going from Ti 2 to V 2 the 
two electrons which are added to the 3d# preferentially couple high spin 
leading to a ground state. When two additional electrons are added 
into the 3d# orbital as in the ground state of O 2 , the overall bonding 

becomes much weaker. This shows that doubly occupying the 3d# orbitals 

* 

Is unfavorable; although part of the difference between V 2 and Cr 2 involves 
3d- 3d exchange interactions. 

One measure of the atomic exchange interaction is ti\e relative energies 
of the atomic states as the coupling of the nd electrons is changed (we 
use Moore[16] to estimate the importance of this effect). For Sc, 4s* 3d 2 , 


the energy needed to change the 3d coupling from 8 F in the 4 F state to 
*D in the 2 D state is 0.8eV. For Ti 4s ! 3d 3 , 4s 1 3d 3 ( 2 G) 8 G is l.OeV above 
4s 1 3d 8 ( 4 F) 5 F. By Cr, with the maximum open 3d orbitals, this value has 
reached 2.5eV. This means that the 3d-3d bonds are expected to be weaker 
for Cr 2 unless the bond length shortens to increase the 3d-3d overlap. 

For the TM elements with more than half filled 3d shells, the formation 
of 3d bonds becomes much less favorable for two reasons. First the ratio of 
4s to 3d orbital sizes is larger for the right half of the first transition ro*. 
Secondly the presence of doubly occupied 3d orbitals leads to repulsive 
interactions which effectively cancel any bonding interactions from the 3d 

* 0 

shell. Thus the bonding here is dominated by the 4s electrons. 

>• 

For the second transition row the qualitative features of the ponding 
from the corresponding three atomic asymptotes are siifiilar. However, the 
more comparable 5s and 4d orbital sizes as well as the smaller 4d exchange 
terms lead to stronger 4d bonding. Hence the 4d£ overlaps appear to 
be more comparable to the 4d o and 4d?r overlaps as evidenced by a low- 
lying state of Nbj with occupation 4dcrj4d?rj4d5j. Also for the right hand 
side of the second transition row the 4d rn + 2 atomic states are available 
for bonding. An additional important difference is that relativistic effects 
become quite important for the right half of the second transition row (see 
Ref. 18, for example). It is for this reason that the current all-electron 

t 

studies are restricted to the lighter elements of the second transition row. 

HL Calculational Details 

As noted in the previous section, both the (n4-l)s 2 nd ro and (n-{~l)s , nd rn "*" 1 
states of the atoms are important in the bonding in the dimers. Hay [19] was 
one of the first to realize that the basis sets optimized for the (n-f-l)s 2 nd m 
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states yield a very poor description of the (n+l)s , nd m + l states, but by op- 
timizing one additional d function for the (n-j-l)s 1 nd m+1 occupation, the 
separation between the (n-f-l)s 2 nd m and (n-j-3)s 1 nd m + 1 occupations was 
very close to that obtained in numerical Hartree-Fock calculations, When 

(n+l)s bonding is present, due to the near degeneracy of the (n-|-l)s and 

• * 

(n-f-l)p orbitals, the basis set requires p functions in the same region, Since 
the (n-f l)p orbital is not occupied, such functions are missing in the atomic 
basis sets. One common approach is to optimize two additional p functions 
for the (n-fljs^n-fljp^d"* occupation and then multiply the exponents of 
the additional functions by from 1.2 to 1.5 to bring them more into the 
region of the (n-f-l)s function [20]. The final concern in the basis set is f 
functions to help describe any d-d bonds. For Cu? with a 4s-4s bond or 
Sc 2 with one-electron d bonds, the addition of f polarization functions 
has very little effect. However, McLean and Liu [21] have shown that 
for Cr 2 , with two-electron d bonding, the addition of f functions lowers - 
the energy by more than leV at the SCF level. Thus it is ‘ important 
to have f functions, especially when comparing one state which has one- 
electron nd bonds with another that has two-electron nd bonds, such as in 
TI2. 


» 

The basis sets used in this work meet these requirements. For the 
first row transition elements the basis set is a (14sllp6d3f)/ [8s6p4d2f] seg- 
mented contraction based on Wachters’ basis set[20] plus his 4p(multiplied 
by 1.5) and the Hsy diffuse 3d {19]. The 4f functions are optimized at the Cl 
level for the 4s 1 3d m + I state of the atom (note that for the Cu 2 calculations 
the three f functions were contracted to one). For the second transition 
row elements the basis set is a (17sl3p9d2f)/[6s5p5dlf] general contraction 
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. based on Huzinaga's basis set[22] plus 5p(mu)tiplied by 1.3), diffuse 4d, and 
4f functions as optimized by Walch, Bauschlicher, and Nelin[23]. The 3s 
function arising from the d, as well as the 4p functions arising from the f 
functions were removed. 

In this work we use the CASSCF approach to perform the MCSCF 

» 

calculations. Ideally one would like to include the nd, (n-j-l)s and (n-)-l)p 
orbitals in the active space to allow the formation of both the s*s bond and 
nd-nd bonds and to account for the (n+l)s to (n-|-l)p near degeneracy. 
More extensive correlation could then be included with a d calculation. 
However, even the CASSCF treatment is not always possible. The first ap- 

* m 

proximation we make is to remove the (n-f-l)p orbital from the active space, 

to 

except for Sc 2 and Y 2 which arise from the (n-f-l)s 2 nd m .^{- (n-|-l)s 1 nd rT, + 1 
atomic limit. Here the (n-J-l)p orbital should be included in the active 
space since the (n-fl)s to (n-|-l)p near degeneracy is important for the 
(n-J-l)s 2 nd r ” state of the atom. For the remaining systems which arise from 
the (n-f l)s 1 nd m + I -f* (n^-l^nd™"*" 1 atomic limit the (n-|-l)p does not 
contribute for the atom and tests show that omitting the (n-fl)p leads to 
less binding energy at the CASSCF level(due to angular correlation effects 
which vanish at large R), but virtually the same D e when correlation is 
added. Thus, this approximation effects only those calculations for which 
only CASSCF calculations were performed. 

The second approximation is to restrict the number of active electrons 
in <r, w and 6 symmetry. The constraints are similar in spirit to those in the 
GVB wavefunction[5]. As an example of the GVB wavefunction consider 
the *£+ state of Ti 2 which has the configuration 

4scrJ 3d o* 3d* 2 . 3dirJ. 
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ie, a quadruple bond with 4so, 3do, and 3d?r bonds. The configurations in 
the GVB wavefunction consist of a gemina) product of the three possible 
configurations constructed from the bonding and antibonding orbitals of 
each pair ie. 81 configurations in the case of four pairs. This wavefunction 

dissociates to neutral atoms and also allows spin recoupling. In the CASSCF 

» 

calculations it is not possible to do precisely the GVB calculation since the 
configuration space must be a product of full Cl calculations within each 
symmetry type. Consistent with these constraints we restrict four electrons 
to the four a orbitals, two electrons to the two n x orbitals, and two electrons 

to the two it, orbitals. The 6 orbitals would not be active in this calculation 

* 0 

since they are forced to have no electrons by our constraints. For the other 

cases the constraints are not specifically stated but can be inferred from 

* , — 

the principle configuration given. 

The last approximation which was made only for the second transition 
row dimers is to freeze the deep core orbitals, the n=l,2,3 shells and 4s 
orbitals at their atomic forms, based on a high spin SCF calculation at 
R=50.0 ao, while the 4p and valence (5s, 5p, and 4d) orbitals are optimized 
at each distance. Since the core orbitals are expected to be atomic like, this 
approximation will have negligible effect on the results, but results in some 
computational savings. 

When it is possible, we follow the CASSCF calculations with multi- 
reference singles and doubles Cl calculations. In these calculations only 
the electrons arising from the (n-f l)s and nd shells of the atom are cor- 
related. The references configurations are selected as the most important 
configurations in the CASSCF wavefunction. In general the most important 
configurations are found to be the leading terms in the GVB gemina] ex- 


* pansion. Different reference lists were used near H, and at long range. The 
long distance points were computed using the high spin SCF wavefunction 
as the reference configuration, while the points near H e use the CASSCF 
wavefunction. 

The calculations for the first transition row dimers use the MOLECULE 
[24]-SWEDEN[25] programs while the calculations for the second transition 
row dimers use BIGGMOLI[26)-SWEDEN. 


Vbs 


IV, First Row TVansition Metal Dimers 

Fig, l shows calculated potential curves for low-lying states of Scg while 
Table D shows calculated spectroscopic constants for all the first row TM 
dimers. The dominant configurations for the states considered here are: 
5 E~ 4scr 2 3daJ 4sai 3d;r‘ y Sd^J* 

8 E~ 4scr% 4s <r* 3d7rJ^ Sdn^ 

4sa 2 4s a\ 3dcrJ 
a E+ 4s<7 2 4s<7 2 3dcrJ 3drri 

The initial study of SC2 by Walch and Bauschlicher[17] found oxJy 
weakly bound states arising out of the 4s 2 3d* -j- 4s 2 3d 1 -asymptote in con- 
trast to mass spectrometric experiments[27] which indicated strong bonding 

r 

(De=l.l±0.2eV), A 8 E^ state was found which was bound by about 0,8 
eV with respect to the 4s 2 3d 1 ■+* 4s 1 4p , 3d i atomic limit, but unbound with 
respect to 4s 2 3d 1 -f 4s 2 3d* . However at about the same time that this work, 
was published, matrix isolation studies by Knight, VanZee, and Weltner[28] 
indicated a bound ®E state of Sc2> "From ESR studies it appeared that this 
state arose from the 4s 2 3d* -f 4s 1 3d 2 atomic asymptote which had not 
been studied in detail in the previous theoretical studies A reexamination 
of this system[30] revealed a new *£7 state which had been missed in the 
previous study because its R e (about 5.0 ao) is much shorter than the R e 
values for the states studied previously (about 7.0 ao). The s E7 state turned 
out to be of considerable theoretical interest because it exhibited multiple 
Sd bonding(three one-electron bonds) and constituted the first theoretical 
evidence of multiple 3d bonding in a first row transition metal dimer. The 
D e given in Table II is obtained from the Cl calculations by computing 
the binding energy of the 5 E7 state with respect to the 4s 2 3d 1 -{-4s 1 3d 2 
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asymptote and then subtracting the experimental asymptotic separation to 
give the D e with respect to 4s 2 3d J + 4s 2 3d 1 . The resulting D e of 0.44 eV 
and cj e of 184 cm -1 are in reasonable accord with experimental values of 
D e =l.leV[27] and u t = 238.9 cm -1 [29). 

Adding two electrons to the 3di$ orbitals of the £ E“ state of Sc 2 and 
high spin coupling the 3d electrons leads to the 7 E+ sUt of Ti 2 . The Ti 2 
7 E+ state, like the SC 2 5 E~ state comes from the 4s 2 3d m -f- 4s 1 3d m + 1 
asymptote. These states have long bond lengths(R e greater than 5.0 &o) 
and small vibrational frequencies^ about 200 cm -1 ) The states of Ti 2 
which are considered here have the configurations: 

7 E+ 4so 2 3doJ 4s a\ Sdn l x% 3 d7rj a 3d^„ 3d5 ( 1 x *_ ya) , 

4s( 7 2 3d(J 2 3d7rJ w 33^ 

* 

For Ti 2 in Table I one sees that the excitation energy to the 4s 1 3d rn + 1 
T-4s 1 3d r7> + 1 atomic asymptote is about half as large as for Sc 2 . Thus- 
the *11+ state which arises from the 4s 1 3d 3 -f- 4s ! 3d 3 atomic asymptote 
becomes a competitor for the ground state of Ti 2 . The bonding here is a 
triple two electron 3d bond(3da, 3d7r 2) and 3d7T y ). This leads to a short 
R state, R e =3.72 ao and u- e =438cm“ 1 . The bond length here is not 
known experimentally but the experimental vibrational frequency is 407.9 
cm" -1 [29] which is consistent with the *£+ state of Ti 2 , hut inconsistent 
with the 7 £+ state. Based on this we tentatively assign the ground state 
of Ti 2 as *£+ although the calculations place *£+ 0.40 eV above 7 £+. 

Table U also shows the calculated D« of Ti 2 . Here we see that CASSCF 
obtains only a small percentage of the binding energy. The Cl wavefunction 
consists of a MRSDCI from nine references (more than 180,000 CSFs). The 
calculated binding energy here is 1.94 eV compared tc *d experimental 



D e of 2,91 eV (31) (all referenced to 4s 1 3d 3 4“ 4s 1 3d 3 ), Thus even when 
more extensive correlation is added, the order of the states is not changed. 
However, Cl calculations describe high spin states better than low spin 
states; the is not as well described as 7 E+. Considering the difficulty 
in describing multiple two-electron 3d bonds, as measured by the errors 
In D e , it would not be surprising to find that a better calculation would 
reverse the order of the states. Thus trends in the correlation treatments 
also support the *E+ assignment of the ground state, 

Fig. 2 shows calculated potential curves for the 3 E“, and l V t states 

of V 2 . The dominant configurations for the states considered here are: 

» • 

3 Ej 4 s< 7^ 3 dOj 3d7r^ w 3d?r^ 3d5^ y ^ 3d(5^3_ y a^ 

3 A, 4s oj 3doJ 3d7r^ 3d7rJ w 3d5^ 3d 6\ ti ^ yi)g *• 

% 

The 3 E“ state of V 2 has the same triple two-electron 3d bond as 
in Ti 2 with the remaining two electrons in the 3d<5 orbitals. Because the* 
3d0 orbitals still have small overlaps in the region near R«(about 3.5 ao), 
the lowest state is a triplet state arising by forming two one-electron 3d5 
bonds. This occupation also gives rise to and *E+ states analogous 
to the corresponding states the wJttJ occupation of 0 2 . The R e and u t 
values obtained from the CASSCF curves[9] are in good agreement with the 
recent results of Langridge-Smith, Morse, Hansen, Smalley, and Merer[32] 
for the 8 E~ ground state of V 2 . The presence of a very low-lying state 
(possibly J E+) was indicated by the spectral analysis of Merer [33] and is 
consistent with the very low-lying l T t state( 1 E+ is obtained as a second 
root of the CASSCF and is also very low-lying). Table H also shows the 
calculated spectroscopic constants for the 8 A # state of V 2 . This state was 
studied because the corresponding state is very low-lying for NT> 2 but here 


i 
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the excitation energy 3 £J~ -► 8 A* is 0.92 eV. 

Fig.3 shows calculated potential curves for O 2 . Here the 3d5 orbitals 
are doubly occupied which is expected to be unfavorable based on the V 2 
result that the 3d <5 orbitals were preferentially singly occupied. This is 
consistent with the weaker bonding in O 2 , D«=1.6eV[34] as compared 
to V 2 , D e =3.0eV[3l]. Because of the weaker bonding in O 2 the CASSCF 
potential curve is not bound|9] .However the potential curve does exhibit a 
shoulder near the experimental R e which is suggestive of an inner well. As 
discussed in section VI a more accurate computational description of the 

quintuple 3d bond in O 2 is a very difficult problem which remains as a 

» • 

major challenge for theory. 

Fig. 4 shows calculated potential curves for TiV'while Table HI sum- 

* » — 

marizes the calculated spectroscopic constants. The 4 H, and 2 A states 
may be derived from the V 2 3 £~ configuration by removing a single electron 
from the 3d<7, 3d;r, and 3dtf orbitals, respectively. The ordering of these ' 
states may be understood on the basis of orbital overlap and intra atomic 
coupling[35] arguments. Based on atomic overlap arguments one expects 
the ordering 

2 A < 4 E“ < 4 H 

However, while all these states dissociate to pure V e D and the 4 S”“ 
and 4 H states dissociate to pure Ti 5 F, the 2 A state dissociates to a mixture 
of 20 % Ti 5 F and 80 % Ti 5 P. It is this latter intra atomic coupling effect 
which leads to the observed ordering of states with 

4 sr < 2 A < 4 H 

The calculated 4 E“ ground state is consistent with the experimental 
results of Van Zee and Weltner which show an ESR spectrum consistent 


• with 4 E""|36). The calculated bond length of TiV is intermediate between 
the bond lengths of Tfe and V 2 . 

Table II also shows results for the *£+ state of CU 2 . The electronic 
configuration here is: 

*E+ 4s 3dcrJ 3daJ 3drrJ 3drrJ 3d<5 4 3d$J 

This configuration arises from the 4s l 3d 10 -f 4s 1 3d 10 asymptote. The 

bonding here appears to involve mainly the 4s electrons with the 3d electrons 

remaining essentially atomic like[37], [38]. However, this viewpoint has been 

challenged by Pauling who argues in favor of 4s 1 4p*3d° hybrids[39]. A 

central concern for the CU 2 calculations has been the calculation of the bond 

length. Here we find that correlation of the 3d electrons shortens the bond 

by 0.19 ao leading to a final bond length of 4.35 ao which is 0.15 ajf longer 

% 

than experiment. It has been argued[37] that the remaining discrepancy is 
largely due to relativistic effects since twice the relativistic contraction * 
of the 4s orbital for the 4s 1 3d 10 state is 0.13 ao. Since that' time cal- 
culations have shown that only about 70% of the atomic core-valence 
contraction in the alkali dimers is observed in the molecular systems[40]. 
Martin[41] using an all electron treatment which includes relativistic effects 
as a perturbation and Laskowsld et.al.[42] who have added relativistic 
effects through an effective core potential approach, both compute the con- 
traction to be 0.10 a c , or about 70% of the atomic contraction. In both 
approaches[41,42] correlation does not effect this conclusion. More recently 
Martin and Werner[43] have considered the effect of higher excitations with 
CEPA, in addition to the relativistic effects, and find a bond length in good 
agreement with experiment. It has been observed[44] that a POLCI treat- 
ment leads to a better D, for CU 2 than a SDCI treatment: However, it is not 
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* clear that such approximations can be used for those systems where there is 
d-d bonding. 
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V. Second Raw TYansition Metal Dimers 

Fig. 5 shows calculated potential curves for selected states of Y 2 . The 
dominant configurations for the states considered here are: 

*E+ 5s cr 2 4d7r 2 u 4d7r 2 u 

5 £J7 5s o\ 5sai 4daj 4d7r* tf 4d7rJ u 

1 A t 5so 2 5so 2 ( 4d7r 2 tl - 4d7r 2 u ) 

*E+ 5s j 2 5so 2 4do 2 

The 5 E“ state here is directly related to the S E7 state of SC 2 and 
arises from the 5s 2 4d* -f- 5s 1 4d 2 atomic asymptote. The long R ! E+ 
and 1 A, states arise from the 5s 2 4d* -f* 5s 2 4d J asymptote. The long R 
*E+ state is directly related to the SC 2 *E+ state, while the *A # state 
was not studied for SC 2 but arises from the same configuration (tt 2 ) as the 
®E“ state of Sc 2 . A new state for Y 2 which was high in energy for Sc 2 
is the short R *E^ state which arises from the 5s 1 4d 2 -f 5s 1 4d 2 atomic, 
asymptote. That this state is low-lying for Y 2 but not for Sc 2 reflects the 
significantly greater 4d bonding in the second transition row. An unusual 
feature of the short R *E+ potential curve is an apparent avoided crossing; 
the state which is crossing the short R ! E+ state is a *E+ state derived 
from the same configuration as the l A g state. This *E+ state should lie 
slightly above the 1 A # curve. An interesting feature of the Y 2 potential 
curves as compared to the Sc 2 potential curves is the more comparable R e 
values for states arising out of the three possible asymptotes. This result 
is consistent with the more comparable (n-f l)s and nd orbital sizes for the 
second transition row as compared to the first transition row. Table IV gives 

i *■ 

calculated spectroscopic constants for the second transition row dimers. No 
experimental data exist for Y 2 . However, Weltner does not observe an ESR 



spectrum for Y 2 as was observed for Sc 2 [45], One possible explanation for 
this is that the 1 £+ state is actually the ground state for Y 2 even though 
the calculations place it 0.87 eV above the statejhere as discussed 
above for Ti 2 we expect more extensive Cl calculations to favor the *£+ 
state over the 5 r~“ state. As was done for the first transition row the T e 
values are obtained by computing the binding energy for each state with 
respect to its corresponding asymptote and then positioning the asymptotes 
at the experimental separations. Thus we are including a correction for 
the error in the Y 2 asymptotic energies due to relativistic effects. Further 
experiments would be necessary to decide definitely whether *£+ is the 
ground state of Y 2 . Here the most useful experimental result would be a 
measurement of H e since R e differs significantly between'the short R *£+ 
and the 5 £~ state whereas the vibrational frequencies are not substantially 
different. 

Fig. 6 shows calculated potential curves for Ntfc while the calculated 
spectroscopic constants are given in. Table IV. The dominant configurations 
for the states considered here are: 

5s(7j 4dtfJ 4d7rJ„ 4d7rJ w 4d$J ff 4dd^ s _ f3 ) # - 
5s a\ 4daJ 4drr^ 4d7rJ u ^\ y9 4d$ ( 1 , 3 _ f3)# 

5saJ 4d a 2 f 4d7rJ 4d 

Here the s £~, l T 9 , *£+, and & A g states are directly related to the 
corresponding states of V 2 . The and 3 4> # states arise by exciting an 
electron from 4d<7 and 4d?r, respectively, into 4d5. From Fig. 6 one sees 

that the electronic structure of seems very similar to that of V 2 in 

» 

that the lowest calculated state is 8 £” and there is a low-lying *r f state. 
There are, however, two additional low-lying states 8 A # *and 8 4> # . Here the 
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. ®A # state is very close in energy to the 3 £“ state while the 3 * f state is 
0.96eV higher, The ordering of these states is consistent with the ordering 
of the 4 £~ and 4 n states of TiV. The increased stability of 3 A f for Nb 2 
as compared to V 2 implies stronger 4d0 bonding in Nb 2 as compared to 
the 3d£ bonding in V 2 . This is also consistent with the stronger bonding 
observed in Nb^ (calculated D,=2.24 eV) as compared to ^calculated 
D«=0.60 eV). 

Experimentally very little is known about the electronic structure of 
Nt> 2 . Smalley has obtained the optical spectrum of but the spectrum 

shows little resolved structure and is too complex to enable analysis[46j. 

* » 

The calculations also show this increased complexity in the Nt 2 spectrum 
due to the presence of a low-lying 3 A, state in addition tcfthe l T t and *£+ 

#, m 9 
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states which were also present for V 2 . Also for Nb relativistic effects are 
much larger than for V and mqy increase the coupling of states. As for 
Y 2 some experimental information would be very useful here especially a 
determination of R € which could help to distinguish between the S A # and 
states as candidates for the ground state of N^. 

In order to compare Cr 2 and Mo 2 the ground *£+ state of M 02 was 
considered at the CASSCF level. This state has the configuration 

5s <rj 4 da* 4dwJ„ 4dwJ. 4 d^ y# 4d^ x9 _ y a)# 
which is directly related to the ground state configuration of Cr 2 . Mo 2 
differs from Cr 2 in that the CASSCF curve is bound due to increased nd 
bonding for the second transition row. The calculated spectroscopic con- 
stants here are(experimental values [47,48] in parenthesis) R e =1.99A(1.93A), 
ct/ e =399cm‘“ 1 (477cm“' 1 ) and D«= 0.77eV(4.2±0.2eV). 

VI. The Description of Electron Correlation in Cr 2 
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As indicated in section IV O 2 might be viewed as a failure of the 
CASSCF procedure since only a shoulder is obtained in the region where one 
expects an inner well(3d bonding region), However, the CASSCF procedure 
with this choice of active space is only a first order description of Cr 2 
which allows the spin recoupling of the 4s and 3d electrons from the large 
R high spin atomic coupling to the region near R e where the system can 
be described formally as a hextuple bond, but does not include all of the 
molecular correlation effects. The particular features of Cr 2 which make 
the description of the region near R e so difficult are i)that we are dealing 

with a quintuple bond and ii) that the 3d overlaps are smaller even near 

* # 

R e than would be the case for say 2p orbitals in a first row molecule. The 

result is that Cr 2 is a molecule which requires a complex MCSCF zero order 

• . • 
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description involving very high orders of excitation withjespect to a single 
configuration SCF reference. 

In order to quantify this somewhat we carried out small valence Cl 
calculations among the active orbitals from the CASSCF calculation(for 
R=3.25ao). Here we allowed single and double excitations from various 
subsets of the most important configurations in the CASSCF wavef unc- 
tion, The first calculation which used as reference configurations the 16 
configurations which are doubles or less from SCF and thus included up to 
quadruples with respect to SCF led to a valence energy 2.42 eV above the 
CASSCF energy. We then added all the quadruples which are important in 
the CASSCF(16 more configurations) so that we are including through hex- 
tuples and obtained an energy 0.70 eV above the CASSCF energy. Finally 
we added the hertuple excitations(6 more configurations) which are impor- 
tant in the CASSCF so that we are now including through octuple excita- 
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• tions and the resulting energy was still 0.42 eY above the CASSCF energy. 
From analysis of the configurations which need to be included to get a 
reasonable valence energy we see that the reference configurations must in- 
clude quadruple excitations among the 3d<5 orbitals and double excitations 
among the 3da and 3d?r orbitals but that the 4s bonding orbital may remain 
inactive(i.e. doubly occupied in all reference configurations). This result 
is in line with the relative overlaps of these orbitals and suggests that the 
3d <5 orbitals are still in the small overlap region even near R e . One also 
expects that V 2 would be easier to describe than Cr 2 since the 3d£ orbitals 
are singly occupied. 

* s 

The CASSCF procedure used does not include all of the important 
molecular correlation effects. For the TM atoms the most important ad- 
ditional correlation effects are derived from the 4p, Sd'fhere 3d' is a tight 
diffuse correlating orbital for the 3d) and 4f orbitals and we expect the most 
important additional molecular correlation terms to involve these orbitals. 
One approach to introducing these missing correlation effects would be the 
direct use of multi configuration singles and doubles CL However, if we 
Hie the CASSCF wavefunction(3088 configurations) as a reference the Cl 
expansion is 57 million configurations for a [8s6p4d2f] basis set. However, 
given the difficulty in reproducing the CASSCF energy for smaller numbers 
of reference configurations it is clear that any Cl calculation with fewer 
references than the CASSCF wavefunction will have to be carried out with 
caution. 

Given these difficulties several alternative methods have been tried for 
introducing the missing molecular correlation in approximate ways. By far 
the most extensive study of additional correlation effects and additivity of 
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correlation effects for Cr 2 is that of Walch[12). In the extended CASSCF 
calculations of Walch, the 3d' shell and 4p shell are added to the CASSCF 
calculation. This approach would be rigorous if one could include all the 
additional active orbitals in one calculation. Unfortunately, this is not pos- 
sible even with orbital occupation constraints and a procedure was used 
in which the extra active orbitals were added by symmetry block and the 
separately obtained energy contributions were summed to obtain an es- 
timate of the extra molecular correlation. The problem with the method 
is that it assumes that the atomic correlation is constant with R, however, 
“atomic correlation” is reduced as the bonds form. For Ti 2 calculations can 
be carried out in which all the 3d' orbitals are included simultaneously and 

it is found that the extended CASSCF procedure overestimates the effect of 

• , ■ 

3d' by a factor of four due mostly to the variation with R of the 3d' atomic 
correlation terms. On the other hand comparisons of extended CASSCF 
calculations and Cl calculations for N 2 and Ti 2 indicate that the approach 
of summing the individual contributions does not give a potential curve 
that is deeper than the Cl curve and that the overestimation of the effect 
of 3d' does not effect the R e or u t significantly. Thus, we conclude that the 
reasonable potential curve obtained for Cr 2 by Walch (see Fig. 3) must arise 
at least in part by cancellation of errors ie. the overestimation of the effect of 
3d' is being compensated for by the omission of other molecular correlation 
effects such as correlation involving 4f. The calculated spectroscopic con- 
stants are (experimental values[34,40-51] in parenthesis) R f =1.78A(1.68A), 
u e =383cm~ 1 (480cm"~ 1 ), and D e =0.71eV(1.56eV). These results, while in 
part due to cancellation of errors, do indicate the importance of 3d' and 4p 
for more extensive Cl calculations. 
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Another approach due to Das and Jaffe|13), is in principle equivalent to 
the extended CASSCF approach used by Walch. These workers use sym- 
metrically orthogonalized atomic orbitals (called PLO’s) in an attempt to 
greatly reduce the size of the MCSCF expansion. Thus, Das and Jaffe did 
not attempt to reproduce the valence correlation in the CASSCF calculations^], 
but rather computed additional correlation effects and then added them to 
our CASSCF potential cum. They predict that 3d' and also 3p shell excita- 
tions contribute to the binding energy. Due to limitations In the Cl expan- 
sion and the small basis used, their conclusions while interesting are some- 
what tentative. 

A computationally very different approach due to Goodgame and Goddard[14] 

r 

assumes that the missing correlation serves mainly to correct the location of 

the ionic atomic asymptotes. Those authors attempt to include these effects 

0 

by empirical modification of the integrals to correct the atomic ionization 
potentials and electron affinities to agree with experiment This method 
does lead to a reasonable potential curve for Cr^ although the bond length is 
somewhat too short which suggests that this method over corrects to some 
extent. 

The main value of these approximate approaches is expected to be to 
lead to a better understanding of the additional correlation effects needed 
to compute a reliable potential cum for Cr 2 . In the opinion of the present 
authors the most promising theoretical method for actually doing the large 
Cl calculations necessary to solve the Cr 2 problem is the externally con- 
tracted Cl method[52]. These calculations will still be difficult due to 
the need for an extensive reference wavefunction leading to a few mil- 
lion configurations. However, it is probable that these calculations will be 
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carried out within the next few years. 

VII. Conclusions. 

The CASSCF procedure is found to generally provide reliable R e and 
w t values. However, D« values are systematically underestimated at the 
CASSCF level; the worst case is Cr 2 where the CASSCF curve is not bound 
although a shoulder is observed in the region near the experimental R<. 
The CASSCF procedure provides a consistent set of calculations for these 
molecules from which qualitative trends and a picture of the electronic 
structure may be derived. 

The CASSCF/CI procedure has lead to confirmation of the 5 E^“ ground 
state of Sc 2 and the ground state of TiV. Numerous predictions remain 
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untested. Among these are the prediction that Tfe has a 4 E+ ground state 
with a short bond length and essentially all the computed results for Y 2 and 
Nb 2 . The computed results for Nl >2 are especially interesting since the ex- . 
perimental spectrum for Nb 2 is so complex that it has not been analyzed. 
Experiment has confirmed the prediction of the ground state symmetry and 
bond length of V 2 . 

In addition to specific predictions, these studies yield a simple qualita- 
tive picture of the bonding in the transition metal dimers. The 4s 2 3d m 
-f* 4s a 3d m asymptote leads to Van der Waals bonding at large R due to 
the 4s->4p near degeneracy effect. Here the 3d electrons are only weakly 
coupled. These weakly bound states have not been observed in experiment, 
and probably are not of chemical interest. If the two atoms have the oc- 
cupation 4s 1 3d m + 1 strong two-electron d-d bonds (ie. .two electrons per 
bond orbital as in Ha, a bond order of one) and a short bond length results. 
Such bonding has been observed for Tfe *£+, Va # £j“, and Cra 


When the excitation energy is very large the cost of promoting both atoms 
to the occupation 4s , 3d m + 1 is so large that bonding arises from the mixed 
asymptote; 4s , 3d m + 1 -f4s a 3d m . Such bonding has been found for Sc 2 C E^ 
and Tia In these states there are three 4s electrons, or a bond order 
of one half. These states also hare long bond lengths, perhaps because of 
the 4s repulsion, which favor one-electron bonds. Such one-electron bonds 
(ie. one electron per bond orbital as in a bond order of one half) and 
are favored at long R. These one-electron bonds naturally retain the atomic 
high spin coupling to maximise the exchange interaction. 

The strengths of the nd bonds are determined largely -by overlap con- 
siderations. Since the relative overlaps of nd orbitals are ndo fz nd?r >ndfi, 
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the ndo and nd?r orbitals are filled before the nd£ orbitals. For the first 
transition row, the d 6 overlap is sufficiently small that two-electron d£ 
bonds do not contribute strongly to the the binding. For example Vj. has . 
both da and d n two-electron bonds and two dS one-electron bonds. In O 2 
where the two additional electrons are added to the d6 orbital, instead of 
increased binding, the D, is smaller for Cr 2 than Va. This arises because 
the small overlap of the d 6 orbitals leads to such weak bonds that they 
do not compensate for the atomic d-d exchange terms which are lost to 
form the bond. The loss of these atomic d-d exchange terms is quite im- 
portant In describing weak bonds. The problems in treating Cr 2 correctly 
are related to this balance of weak two electron d 6 bonds and the atomic 
d-d exchange interactions. Since this problem is not present in the other 
dimers considered, our CASSCF treatment fails only for O 2 . 

There are several differences between the first and second transition 
row: i)<r*«>/<r 4 rf> is smaller than <tu>/<tm> for comparable loca- 


lions in the first row ii) the exchange interactions are smaller for 4d than 
for 3d and iii) the 5s 2 4d m -^5s 1 4d m + 1 excitation energy is smaller than the 
4s 2 3d m -♦45 , 3d ,r, + , excitation energy. These features lead to stronger 4d 
bonding than 3d, with a notable increase in the strength of the d5 bond and 
a reduced significance of the mixed asymptote. There are several examples 
of these changes. Ya has a low-lying *£+ state, which is a very likely 
candidate for the ground state, arising from 5s 1 4d 2 + 5s 1 4d 2 . M 02 has a 
larger D, than O 2 and the CA&SCF treatment is now able to describe this 
bonding. Nba has a low-lying 8 A, state (4dcrj4dirj4d5j) which arises by 
moving a 4d*r electron into a 4d£ orbital, while the corresponding state is 
0.92 eV up for V 2 
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Tabla I a. Ralativa orbital alias and azcltatlon anarglaa lor fir at row 
transition matal atoms. Tha orbital azpactatlon valuta art takan 
from Numarlcal Hartraa-Fock calculation, whila tha azcltatlon 
anarglas art takan from azparlmant. 


Ato^ <r 4 #>/<rsa> 4s 2 3d m -*4a , 3d m + 1 4a 2 3d m ->3d fn + 2 


Sc 

2.03 

1.43 

4.10 

T1 

2.32 

0.61 

3.35 

V 

2.51 

0.25 

2.47 

Cr 

2.60 

-1.00 

3.40 

Mn 

2.81 

2.14 

5.50 

Fa 

2.05 

0.87 

4.07 

Co 

3.06 

0.42 

3.36 

11 

3.22 

-0.03 

1.71 

Cu 

3.36 

-1.40 

>• 
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Tabls lb. RtlatlT* orbital slzss and sxcitatlon snsrgiss for sscond row 
transition astal atoms . Tbs orbital szpsctatlon valuos art taksn 
from Nuasrical Hartrss-Fock calculation, vblls tbs sxcitatlon 
snsrglss ars taktn from oxpsrinsnt. 


Atom <rj, >/<Tu> 5« 2 4d m -.5* 1 4d m + I 5« 2 4d ra -.4d m + 2 


Y 

1.61 

1.36 

3.63 

Zr 

1.79 

0.59 

2.66 

Mb 

1.92 

-0.18 

1.14 

Mo 

2.05 

-1.47 

1.71 

Tc 

2.16 

0.41 


Ru 

2.29 

-0.87 

0.22 

Rb 

2.42 

-1.63 

-1.29 

Pd 

2.54 

-2.43 

-3.38 

Ag 

2.67 

-3.97 

>• 
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Tabla II. Calculatad spactroscoplc constants lor first rov 
transition matal dinars 


Molacula 

Re (A) 

w t (cm" 1 ) 

D c (aV) 

T e (aV) 

Bc 2 5 Er 

Cl 

Exp. 

2.70 

184. 

238.0° 

0.44 
1 .13* 


Il2 l Z+ 
CASSCF 
Cl 
Exp. 

1.07 

1.07 

436. 

438. 

407.0° 

0.37 

1.04 

2.01 c 

0 0 

0.40 

Il2 7 S+ 
Cl 

2.63 

205. 

x.st 
• , 

0.00 

A 

V 2 S E- 

CASSCF 

Exp. 

1.76 

1.76 rf 

564. 

537.5° 

% 

0.60 ' 
2.29*.2.97 c 

0.00 

v 2 ‘ r , 

CASSCF 

1.80 

486. 

0.58 

0.02 

V 2 3 A, 
CASSCF 

1.77 

413. 

-0.32 

0.02 

Cr 2 »£+ 
CASBCF 
Exp . 

1.68 e 

480.° 

-1.4 

1.56^ 



* Raf. 20 . 

1 Raf. 2T 

* Raf. 31 
4 Raf. 32 

* Raf. 49-51 
f Raf. 34 
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Table HI. Calculated spectroscopic constants for the T1V nolecule 


State 

Re (A) 

U e (C *~ 1 ) 

D e (eV) 

4 ir 

1.66 

405. 

0.80 

4 n 

1.05 

410. 

0.15 

2 a 

1.80 

472. 

0.43 


Calculated with respect to TlUs’Sd 3 ) + VCds’dd 4 ). 
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T&bli IV. Calculittd iptctroicoplc con* tint* for Bicond 
transition Bttil dlatrs 

Molicult 

Rf (A) 

U t (CB” 1 ) 

D e <iV) 

y 2 5 e~ 

CI 

3.03 

171. 

2.44 

r 2 1 E+ 

CASBCF 

CI 

2.73 

2.74 

205. 

206. 

1.74 

2.93 

Kba 8 E f - 

CAS8CF 

2.10 

446. 

2.24. . 

m >2 ‘r, 

CASBCF 

2.11 

427. 

2.1V 

»3 *A ( 
CASSCF 

2.01 

501. 

% 

2.12 

n>2 3 *, 

CASSCF 

2.19 

340. 

1.26 

Moa >£+ 
CASSCF 

Exp. 

1.99 

1.93° 

399. 

477.® 

0.77 

4.2* 


• Bif. 47 
1 Rif. 46 


3C 


row 
T< (#V) 
0.00 

0.87 

0.00 

.0.09 

0.12 

0.96 
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Figure Captions. 

Fig. 1, Calculated CASSCF/CI potential curves for selected states of Sc 2 . 
Fig. 2. Calculated CASSCF potential curves for selected states of V 2 . 

Fig. 3. Calculated potential curves for the *£+ state of Cr 2 - The curve 
labeled CASSCF + 3d' includes an estimate of 3d' correlation from the 
extended CASSCF calculations of Walch (Ref. 12);while the curve labeled 
CASSCF -f- 3d' 4 * 4p(Est.) also includes an estimate of the 4p correlation 
based on extended CASSCF calculations on Ti 2 (Ref 12.). Because the 
correlation contributions are obtained by summing individual energy con- 
tributions the effect of 3d' is expected to be significantly overestimated (see 
text). Thus, we expect that a CASSCF calculation in which all the 3d' 
orbitals were added simultaneously would show a much reduced binding 
energy contribution due to 3d' leading to a potential curve with a more 
distinct shoulder than the CASSCF curve but no well. 

Fig. 4. Calculated CASSCF potential curves for selected states of TiV. 
Fig. 5. Calculated CASSCF/CI potential curves for selected states of Y 2 . 
Fig. 6. Calculated CASSCF potential curves for selected states of Nb 2 . 
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